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qq ' The exchange interaction of polaronic carriers with localized spins leads to a ferromag- 

, netic/paramagnetic transition in doped charge-transfer insulators with strong electron- phonon cou- 

^\ • pling. The relative strength of the exchange and electron-phonon interactions determines whether 

the transition is first or second order. A giant drop in the number of current carriers during the 
q , transition, which is a consequence of local bound pair (bipolaron) formation in the paramagnetic 

phase, is extremely sensitive to an external magnetic field. Below the critical temperature of the 
transition, T c , the binding of the polarons into immobile pairs competes with the ferromagnetic 
exchange between polarons and the localized spins on Mn ions, which tends to align the polaron 
moments and, therefore, breaks up those pairs. The number of carriers abruptly increases below 
T c leading to a sudden drop in resistivity. We show that the carrier density collapse describes the 
colossal magnetoresistance of doped manganites close to the transition. Below T c , transport occurs 
by polaronic tunneling, whereas at high temperatures the transport is by hopping processes. The 
transition is accompanied by a spike in the specific heat, as experimentally observed. The gap fea- 
ture in tunneling spectroscopy is related to the bipolaron binding energy, which depends on the ion 
mass. This dependence explains the giant isotope effect of the magnetization and resistivity upon 
substitution of ls O by ls O. It is shown also that the localization of polaronic carriers by disorder 
cannot explain the observed huge sensitivity of the transport properties to the magnetic field in 
doped manganites. 
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I. INTRODUCTION 

The existence of a metal-insulator transition in lanthanum manganites was established in the early 1950s [Q and has 
ly-^ ■ been extensively studied thereafter. The transition is associated with unusual transport properties, including large 
£>T) [ magnetoresistance in the vicinity of the transition, studied in a family of doped manganites with perovskite structure 
£N| ■ with the chemical formula Rei-^D^MnOa, where Re is the rare earth (Re = La, Pr, Nd), and D is the divalent metal 
(D = Ca, Sr, Ba). It is worth mentioning the early studies of the transition in Lai^Pb^MnOa J2], followed by the 
OO ! studies of Prx^Ca^MnOa §, Nd .5Pb . 5 MnO3 §, La . 67 Ba . 3 3MnO 3 g, Lao. 75 Cao. 25 Mn03 ff, La 1 _ I Ca 1 Mn0 3 
(see review ||). The recent resurgence of interest in these systems is related to the demonstration of a very 
large negative magnetoresistance in thin films |5|,[7j [sometimes termed colossal magnetoresistance (CMR)], which 
immediately raised the possibility of technological applications. The colossal magnetoresistance is not limited to 
^ doped perovskite manganites, but was also observed in pyrochlore manganites, chromium spinels and some other 
■ systems, like europium compounds. 

The metal-insulator transition in lanthanum manganites Jl],0,D has been traditionally attributed to a 'double 
O ■ exchange' mechanism, which results in a varying band width of holes doped into the Mn 3+ (i-shell as a function 
J-^ of the doping concentration and temperature JhJ. Recently it has been realized jy], however, that the effective 
carrier-spin exchange interaction of the double-exchange model is too weak to lead to a significant reduction of the 
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electron bandwidth and, therefore, cannot account for the observed scattering rate 12 (see also Ref. 113] ) or for 



localization induced by slowly fluctuating spin configurations |14| . In view of this severe shortcoming of the double 
■ exchange model, it has been suggested [|ll| that the essential physics of perovskite manganites lies in the strong 
coupling of carriers to the Jahn- Teller lattice distortion. The argument [[tlj was that in the high-temperature state 
the electron-phonon coupling constant A is large (so that the carriers are polarons); as temperature decreases the 
growing ferromagnetic order increases the bandwidth and thus decreases A sufficiently for metallic behavior to occur 
below the Curie temperature T Cl in accordance with polaron theory |0|. A giant isotope effect the sign anomaly 
of the Hall effect, and the Arrhenius behavior of the drift and Hall mobilities |l7]] over a temperature range from 
2T e to 4T C unambiguously confirmed the polaronic nature of the carriers in manganites. Polaron hopping transport 
accounts satisfactorily for the resistivity in the paramagnetic phase JlTf] , 

However, the known relation between magnetization and transport below T c and the unusual magnetic ion dynamics 
have prompted the conclusion that polaronic hopping is also the prevalent conduction mechanism below T c urn . Low- 
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temperature optical |19|-|21[], electron-energy- loss (EELS) J22| and photoemission spectroscopies |23j showed that the 
idea [fyjjlj] of a 'metalization' of manganites below T c is not tenable. A broad incoherent spectral feature 19-21 23f| 
and a pseudogap in the excitation spectrum 23 -p5 were observed while the coherent Drude weight appeared to be 
two orders of magnitude smaller |]2Cf| than is expected for a metal, or even zero in the case of layered manganites 
[ p3| . EELS [^2| confirmed that manganites are charge-transfer doped insulators having p- holes as the current carriers 
rather than d (Mn 3+ ) electrons. The photoemission and O Is x-ray absorption spectroscopy of Lai-^Sr^MnOa showed 
that the itinerant holes doped into LaMn03 are indeed of oxygen p character, and their coupling with the d A local 
moments on Mn 3+ ions aligns the moments ferromagnetically p^| . Moreover, measurements of the mobilit y |9|j27[| 
do not show any field dependence and there are significant deviations from Arrhenius behavior close to T c p8| , |17| . 
The resistivity calculated from the modified double-exchange theory is in poor agreement with the data and the 
characteristic theoretical field (~15T) for CMR is too high compared with the experimental one (^4T) jll]]. As a 
result, self-trapping above T c and the idea of metalization below T c do not explain CMR either. Carriers retain their 
polaronic character well below T c , as manifested also in the measurements of resistivity and thermoelectric power 
under pressure Q • 

Therefore, the experimental evidence overwhelmingly suggests that the low-temperature phase of the doped man- 
ganites is not a metal, but a doped polaronic semiconductor. The double exchange and the presence of polaronic 
carriers are insufficient to explain the physics of colossal magnetoresistance. One can also add that there are known 
classes of CMR materials where it is guaranteed that double exchange is non-existent, like in pyrochlore manganites, 
chromium spinels and other compounds. 

In the present paper, we propose a new theory of the ferromagnetic/paramagnetic phase transition accompanied 
by a current carrier density collapse (CCDC) and CMR. Taking into account the tendency of polarons to form local 
bound pairs (bipolarons) as well as the exchange interaction of p polaronic holes with d electrons, we find a novel 
ferromagnetic transition driven by non-degenerate polarons in doped charge-transfer magnetic insulators. The crux of 
the matter is that in the paramagnetic state above the critical temperature a large fraction of polarons is bound into 
immobile pairs (bipolarons). As the temperature decreases in the paramagnetic phase (T > T c ), so does the density 
of mobile polarons, and the resistivity quickly increases with the decline of the number of carriers. With the onset of 
ferromagnetic order at T c , the situation changes dramatically. As a result of the exchange interaction with the localized 
Mn spins, the energy of one of the polaron spin sub-bands sinks abruptly below the energy of the bound pairs. The 
pairs break up, the density of carriers (mobile polarons) jumps up, and the resistivity suddenly declines, as observed 
experimentally. The occurrence of the deep minimum in the carrier density close to the transition point, which 
we suggest calling a current carrier density collapse, allows us to explain the magnetization and temperature/field 
dependence of the resistivity of Lai-xCa^MnOa close to T c as well as the giant isotope effect, the unusual tunneling 
gap, and the specific heat anomaly. 

II. FERROMAGNETIC TRANSITION IN DOPED MANGANITES 

The Hamiltonian containing the physics compatible with the experimental observations mentioned above is 

k,s k,j 
k,q,s 

+ 5> Wq (& q 6 q + l/2), (1) 
q 

where is the local density approximation(LDA) energy dispersion |^0| , h^ s is the annihilation hole operator of a 
(degenerate) p oxygen band with spins s =f and J,, J p( i is the exchange interaction of p holes with four d electrons of 
the Mn 3+ ion at the site j, = h^ h^ — h^h^i, S* is the z-component of Mn 3+ spin, which is S — 2 due to the 
strong Hund coupling, 7iHund, of the four (i-electrons on Mn 3+ sites, N is the number of unit cells. The two last terms 
of the Hamiltonian describe the coupling of p holes with phonons and the phonon energy, respectively (7 q = ~7l q 

is the coupling constant [^5|). The Hamiltonian also contains spin-flip processes, H s f, like h\.,^h^ + h.c, and 

terms with non-diagonal components of the polaron magnetization operator mfk = h^^h^ — h^Ji^i, which are not 
essential for our discussion. If the holes were doped into the d shell instead of the p shell, the Hamiltonian would be 
similar to the Holstein t — J model f3l| with about the same physics of CMR as proposed below. 
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The essential results are readily obtained within the Hartree-Fock approach for the exchange interaction |32|] and 
the Lang-Firsov polaron transformation []33f which removes terms of first order in the electron-phonon interaction in 
Eq. (|), H = e u He- u , where 

U = J2 h lh ls u m (bl cl + b (l ), (2) 

hj S = N^ 1 / 2 ^ k hk s exp(zk • Rj), ujq = (2iV) -1 / 2 7 q exp(iq • R^-), and R^ is the lattice vector. 
With the use of this transformation one finds spin-polarized p bands 

e kT(l) = {+)\j P dS(T - {+)^ B H. (3) 

Here 

£ k = Jj E Use^-^e-ti « E^e-v 2 (4) 

where 



E n 2 1 1 - cosf i • ( R * - r j)] coth 

q 



= ™ E N 2 f 1 - cos 1 ' ( R * - R J')1 coth ( ^ J (5) 



where g 2 ~ 7 2 is the characteristic value of a?- and w q is the phonon frequency. Equation (||) describes the polaronic 
band narrowing |53| and the isotope effect JL6[ . The bare hopping integrals tu define the unrenormalized LDA (local 
density approximation) band dispersion in the initial Hamiltonian ([!]) E^ = Y^i j hj exp[ik • (Rj — Rj)]- In Eq. (||) 
er is the normalized thermal average of the Mn spin, found from the equations below; H is the external magnetic field, 
and fis is the Bohr magneton. The p — d exchange interaction depends only on total (average) magnetization because 
we assume that the system is homogeneous. In addition to band narrowing, the bands shift rigidly down, the value 
of this polaron shift E p is 

£p = 2^E^M 2 - (6) 

q 

The ions Mn 3+ are subject to a molecular field Jpdm/ (2gMn^B), according to (Q), and their magnetization a = 
(S*)/S is given by 

D { J P dm + 2g Mn n B H\ 
a = B s[ 7^-7? . (7) 



2k B T J 
with m the magnetization of holes determined as 

m= 1 £<m k ) = / deN^(e) [/ p (e kT ) - f p (e H )] . (8) 

k J 

Here B s (x) = [1 + 1/(25)1 coth [(5 + l/2)x] - [1/(25)] coth(ar/2) is the Brillouin function, g Mn the Lande g-factor for 
Mn 3+ in a manganite, N( p '(e) the density of states in the narrow polaron band, and / p (eks) = exp(ei cs /fcBT') + l] _1 
the Fermi-Dirac distribution function with y = cxpip/ksT) determined by the chemical potential [i. Note that for 
Jpd < (antiferromagnetic coupling) the main system of equations (|l6|)-(|l8|) remains the same after a substitution 

Jpd I Jpd | ■ 

Along with the band narrowing effect, the strong electron-phonon interaction binds two polarons into a local pair 
(bipolaron), as described in detail in Ref. Jl5| . These bipolarons are practically immobile in manganites because of 
the strong electron-phonon interaction, in contrast with cuprates, where bipolarons are mobile and responsible for 
in-plane transport Q , owing to their geometry Q and their moderate coupling with phonons Q . 

If these bound pairs are extremely local objects, i.e. two holes on the same oxygen, then they will form a singlet. 
If, however, these holes are localized on different oxygens, then they may well have parallel spins and form a triplet 
state. The latter is separated from the singlet state by some exchange energy J st , with some interesting consequences 
discussed below. Because of their zero spin, the only role of the singlet bipolarons in manganites is to determine the 
chemical potential fj,, which can be found with the use of the total doping density per cell x [p4[. 
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The interplay between the localization of p-holes into bipolaron pairs and the exchange interaction with the Mn d A 
local moments is responsible for CMR. The density of these pairs has a sharp peak at the ferromagnetic transition when 
the system is cooled down through the critical temperature T c . As the system is cooled, but is still in a paramagnetic 
state above T c , an increasing fraction of the polarons forms immobile pairs (bipolarons), and the resistivity of the 
system increases. Below T c , the binding of polarons into immobile pairs competes with the ferromagnetic exchange, 
which tends to align the polaron moments and, therefore, breaks those pairs apart. The number of carriers abruptly 
increases below T c leading to a sudden drop in resistivity. These competing interactions lead to the unusual behavior 
of CMR materials and the extreme sensitivity of their transport to external fields. 

To prove the point, we shall find the thermodynamic potential and solve for its extremal value to find the equation 
of state for the polarons. The thermodynamic potential f2 

= fi p + Obp + &s + ^JpdSam, (9) 

has contributions from polarons, bipolarons, localized Mn 3+ spins, and the double-counting term, respectively. For 
the polarons 

fl p = -k B T J deN^ (e) In (l + ye- £ / fcsT ) , (10) 

s 

where iV« (e) is the density of spin polarized states in the polaron band, and y = exp(n/k B T). We can easily 
estimate the integral, assuming that the critical temperature of the ferromagnetic transition is comparable with the 
polaron and bipolaron bandwidth |37|. Then (bi)polarons are not degenerate in the relevant temperature range, 

(A — e)/k B T], and we get 



/ p ~ yexp(-e/k B T) and / bp ~ y 2 exp 



» p = -2^rcosh J ^+f Bg , (11) 



where i>(— 3) is the degeneracy of the polaron p band. 

Polarons, bound in bipolarons with a binding energy A, give a contribution 



ftbp = -k B T In (l + V 2 y 2 De^ kBT ) , (12) 



where D accounts for the presence of triplet bipolarons (see below Sec. VI). We shall consider here a simple case when 
the separation of the triplets from the singlets, J S i, is much larger than the critical temperature. In this case D = 1. 
Finally, for the localized spin contribution we will have 

sinh(5+ \)ri , 

{l s = -k B T\n } 2 " , (13) 

sinn 

with X] = {\J pd m + g Mn ^BH)/k B T. 

The density of polarons n = — {d£l v /d[i)T is found from the condition that the total number of carriers is given by 
the doping concentration x p4| : 

x = -{dn/dn) T , (14) 

whereas one can find equations for the magnetization and the normalized spin a from the following conditions: 

(dn/d<r) T = (dn/dm) T = 0. (15) 

Thus, we obtain the following main system of mean field equations, assuming for a moment that the contribution 
from triplet bipolarons is small (D = 1): 

n = 2vycosh[(a + h)/t], (16) 
m = ntanh[(<7 + h)/t], (17) 
a = B 2 [(m + ih)/(2t)], (18) 

and 
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FIG. 1. Schematic of free polaron (P) and polaron bound pair (BP) densities of states at temperatures below and above T c 
for up (|) and down (J.) spin moments. The pairs (BP) break below T c if the exchange J v dS between p-hole polarons and Mn 
d 4 local spins exceeds the pair binding energy A, as in the case shown. The exchange interaction of polarons with the localized 
spins sets in below T c , the spin-up polaron sub-band sinks abruptly below the bipolaron band, causing the break-up of the 
immobile bipolarons (left panel). A sudden drop (collapse) of the density of the current carriers (polarons) in the vicinity of 
the ferromagnetic transition is the cause of a large peak in resistivity and colossal magnetoresistance. 



= cxp (-26/t), (19) 

which follows from ([lJJ). When triplet bipolarons become important, one should replace Eq. ([l9]) by the more 
accurate Eq. (p8|). Here we use the dimensionless temperature t = 2ksT '/ '(J p dS), magnetic field h = 2/j,bH/(J pi iS), 
the bipolaron binding energy 5 = A/(J p dS), and v{= 3) is defined after Eq. (fTl|). 

The polaron density n is determined by Eqs. (|l6| ), ( |Ts| ) with a — above T c . At the critical temperature, the 
polaron density has a minimal value n c ~ (2x) 1 ^^exp(—S/t c ), it then grows with temperature and saturates at 
n = (1 + 2a;) 1 / 2 — 1 at large temperatures. This is reminiscent of ordinary semiconductor behavior. 



III. COMPETING INTERACTIONS AND CARRIER DENSITY COLLAPSE 



The remarkable observation is that there is a sharp increase of the polaron density (and the conductivity) at 
temperatures below T c . The polaron density approaches the total density x at T — ► if S = A/J p dS < 1, as one can 
see from Eq. (|l6| ) with a saturated magnetization a = 1. The physical origin of the unusual minimum of the current 
carrier density at T c lies in the instability of bipolarons below T c due to the exchange interaction of polarons with 
Mn d electrons. The spin-polarized polaron band falls below the bipolaron band, so that all carriers are unpaired at 
T = if JpdS > A. The evolution of the Hartree-Fock bands with temperature, which corresponds to this behavior, 
is illustrated in Fig. [l[ Note that at all T > T c the position of the polaron bands is fixed at A/2 above the bipolaron 
band, since there a = rn = (|[). Their population depends on temperature via the chemical potential. The exchange 
interaction of polarons with the localized spins sets in at T c , and in the low-temperature ferromagnetic phase one of 
the polaron spin sub-bands sinks abruptly below the bipolaron band, causing the break-up of the immobile bipolarons. 
This interesting feedback mechanism can result in either a continuous or discontinuous ferromagnetic transition, as 
follows from a simple analysis below. 

Linearizing Eqs. (fl6|)-(|l8|) with respect to a and m near T c , we find the critical temperature in zero magnetic field 

t c = K/2) 1 / 2 , (20) 
where the polaron density at the transition n c is determined by 
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FIG. 2. Inverse polaron density x/n in a doped charge-transfer insulator for different magnetic fields h = q^lbH / J p dS , 
A/JpdS = 0.5, doping x = 0.25. A is the pair binding energy, J p dS is the exchange energy of the O p hole polarons with Mn d 
localized spins. For other notations see text. Note that the transition is a strong first order, and then becomes continuous when 
the external magnetic field exceeds some critical value. Inset: temperature of the phase transition as a function of external 
magnetic field. 



n l/2 ln 2 ( X O =2 3/3ft (21) 

K 

It is easy to see that this transcendental equation has solutions only for S below some critical value 5 c (x). This 
means that for 6 > S c (x) the ferromagnetic phase transition is first order with jumps of the polaron density and the 
magnetization (38), as has been observed |^| . The transition is continuous when <5 < S c (x). The numerical solution 
of the system Eqs. (|l^)-((l^) defines the crossover between first- and second-order phase transitions |S8| . 

A relatively weak magnetic field has a drastic effect on the inverse carrier density, 1/n, near the first order phase 
transition, or second order phase transition close to first order. A field equal to only 0.005 J p dS / (2 ps) reduces the 
carrier density collapse by more than a factor of two, and a field of 0.01J pc iS / '(2/jb) changes the transition into the 
continuous one, Fig. ^|. This behavior directly relates to the colossal magnetoresistance found in doped manganites, 
as we shall discuss in the following section. 

One can draw an analogy of this situation with singlet magnetism, e.g. in Pr compounds ]4(i| ] . In this case the ground 
state of magnetic ions is singlet. Depending on the ratio between the exchange constant and singlet-triplet (doublet) 
energy gap produced by crystal-field splitting, there exist first- or second-order phase transitions into a ferromagnetic 
state. In our case the triplet states become important when J st < A, the larger statistical weight of triplet bipolarons 
leads to a deeper minimum in the density of polarons at the critical temperature, and, therefore, to a larger jump in 
resistivity. The effect of polaron binding in a triplet state will be discussed below. 



IV. COLOSSAL MAGNETORESISTANCE 



As a result of the carrier density collapse, the resistivity p — l/(ert/i p ) has a sharp maximum, which is extremely 
sensitive to the magnetic field in the vicinity of T c . In fact, our theory, Eqs. (|l^)-([l8|), describes all the major features 
of the temperature/field dependence of p{T) with a temperature and the field independent polaron drift mobility 
[i p in the experimental range of the magnetic field, Fig. p(a),(b). It gives the correct magnitude of the effect on 
resistivity and explains the extreme sensitivity to external magnetic fields. This suggests that current carrier density 
collapse is the origin of CMR. 

In general, one has to take into account the temperature dependence of the polaron mobility to extend our theory 
for temperatures far away from the transition. At high temperatures, the mobility fi p of polarons is dominated by 
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FIG. 3. Resistivity of Lao.7sCao.2sMn03 calculated within the present theory for A = 900K, J v dS = 2250K for a temperature 
independent mobility (a). The experimental results [8] are shown on panel (b). Note the extreme sensitivity of the theoretical 
resistivity to the external magnetic field (a) , also observed experimentally for the doped manganite (b) (thin solid line is a guide 
to the eye). Panel (c): Resistivity calculated with a temperature dependent mobility according to Eq. (26) with loq = 50meV 
and E a = 300meV and temperature dependent polaron density from panel (a) compared with the experimental results. Note 
the crossover of the transport mechanism from low-temperature tunneling to high-temperature hopping at about the transition 
temperature. For notations see caption to Fig. El 



hopping events since the polaron narrowing factor g 2 grows linearly with T, making tunneling in a narrow polaron 
band virtually impossible at fcgT > hwo/2, where ujo is the characteristic phonon frequency |33| . A simple estimate 
for the so-called adiabatic hopping conductivity together with the Einstein relation between diffusion constant and 
mobility immediately yields 

^-^^cM-E a /k B T), (22) 

where /io = ea 2 /% is the characteristic mobility (one can estimate a as the 0-0 distance in manganites), and E a is 
the activation energy for the hopping. Tunneling mobility is given by 

< n) = "«& < 23 > 

with the relaxation time r estimated by Lang and Firsov f33fl 

r « (E a /t) 4 {Au;/u, 2 } smh 2 (hio /2k B T), (24) 
where t is the characteristic bare hopping integral tij , and Alo is the phonon dispersion. The resistivity is then given 

by 



P 



= 1/ CT , a = ne(M (tun) + M (h ° p) )- (25) 



With our low polaron density at the transition the polaron mobility is = 0.2 cm 2 /Vs for x — 0.33 |l7| , and about 
0.03 cm 2 /Vs for x = 0.25 M that lies in the range typical of polaronic conductors like Ti02 at room temperature 
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p5[ . We have fitted the observed resistivity to the above expression [Fig. g(c)] using for uiq a value of 50 meV, which 
is close to the phonon cutoff in LCMO (50-70 meV @). 

The fit indicates that the activation energy is close to E a = 300 meV. A crossover from tunneling to hopping occurs 
around the critical temperature T c , which is not very different from huo/2kB [S3]. Agreement with the experiment 
(Fig. |^(c)) supports the idea that the temperature dependence of the resistivity is due primarily to CCDC. The 
temperature dependence of the small polaron mobility then allows the resistivity far away from the transition both 
above and below T c to be explained. 



V. ANOMALOUS SPECIFIC HEAT 



The carrier density collapse is also evident through anomalies in thermodynamic quantities. Indeed, we have 
shown above that the ferromagnetic transition is first order, or second order close to first order, as observed. The 
thermodynamic potential changes rather abruptly in the vicinity of the phase transition and this results in a sharp 
peak in the specific heat C, Fig. [|, which has been observed [Q. Note that this is not a result of critical fluctuations 
as suggested earlier |j43| , since they are absent or severely suppressed when the phase transition is first order, or close 
to it. We see that our theory is in quantitative agreement with the experiment for this anomalous thermodynamic 
quantity. 
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FIG. 4. Calculated anomalous part of the specific heat for different values of the magnetic field H. Inset: experimental 
results for Lao.67Cao.33Mn03 [42] (thin solid line is a guide to the eye). 



VI. TRIPLET BIPOLARONS 



Let us now discuss the modification which arises if we include exchange between O-holes bound into bipolarons. 
This exchange generally induces a splitting J s t between singlet and triplet states of the bipolaron. This changes 
somewhat the thermodynamic potential of the bipolarons, since the triplet is subject to a Zeeman splitting. The 
factor D is then 

D = l + e- J * t/fcBT sinh(3£/2)/sinh(£/2) (26) 

as it accounts for thermal excitations of singlet bipolarons into the triplet state, separated from the singlet by the 
energy J st . The parameter 

£ = (J pd Sa + V bp m + gn B H)/k B T, (27) 
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FIG. 5. Inverse polaron density x/n for different magnetic fields for a system with triplet and singlet bipolarons versus 
temperature t = 2ksT / J p dS ■ (A/J p dS = 0.5, doping x = 0.25, and we assume J 3t <C A). The jump in carrier density is 
much larger in a system with triplet bipolarons, but the critical temperature and sensitivity of the critical temperature to the 
magnetic field is lower in comparison with singlet bipolarons. For notations see caption to Fig. 



depends on exchange interaction of the bipolarons with Mn 3+ spins given by the exchange constant J p d, and delocalized 
polarons, given by the exchange constant Vb p . Note that D = 4 at J^/Zc^T <C 1, whereas D — 1 for J st /ksT ^> 1, 
which reflects the higher statistical weight of triplet states compared to singlets. 

It is assumed, as is usually the case, that the triplet states lie higher in energy than the singlet state, J s t > 0. If 
the singlet-triplet splitting becomes smaller than the gap, J s t ^ A, then, because of a higher number of the triplet 
states, their thermal population leads to a deeper minimum in the density of polarons and, therefore, to a larger 
jump in resistivity (Fig. |5|). The dependence of the population of the triplet states on external field makes the system 
somewhat less sensitive to the field. We make an essential assumption that J s t > and that the exchange between 
spins on Mn and triplet bipolarons, J p d, is suppressed to values <C J p d because the bipolarons are strongly localized 
[we also expect that the exchange constant Vb p is the smallest one in (p7|)]. Otherwise, the triplet bound pairs, if 
they were formed in the paramagnetic phase, can survive in the ferromagnetic phase thus reducing or eliminating the 
carrier density collapse. 

The equation (|l^) is changed to read 

whereas the main system of equations (p^|)-(pj|) remains the same. D is given by Eq. (^6|). The effect of triplet 
bipolarons on the thermodynamics of doped manganites becomes insignificant when J s t/A > 1, and the results are 
similar to the case when only singlet bipolarons are involved. The polaron density at the transition n c is determined 

by 

n l J 2 In - '^"^ m = 2*/% (29) 
Dn£[l - (x - n c )/2\ 

which is similar to the case of singlet polarons and also indicates a crossover from first- to second-order phase transition. 
We compare the carrier density collapse in a system with triplet bipolarons to that with singlet bipolarons alone in 
Fig. |^. The jump in the carrier density at the transition is a few times larger in this case as compared to singlet 
bipolarons. At the same time the critical temperature shifts to lower values, and the sensitivity to external magnetic 
field slightly reduces. 
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FIG. 6. Relation between the gap 8 = A/J p dS and the critical temperature t c = 2fcsT c / J v dS calculated from the present 
theory. Inset: tunneling gap in the density of states for samples with different temperatures of the transition: Lao.sCao.2MnC>3, 
2t r =196K; (NdLa) .73Pbo.27Mn0 3 (Tt r =275K); Lao. 7 Pbo.3Mn0 3 (T tr =338K) [24]. For notations see caption to Fig. [| 



VII. TUNNELING GAP AND GIANT ISOTOPE EFFECT 



Recent tunneling measurements have shown that in the vicinity of T c a gap in the quasiparticle spectrum open s up 
[p4|]25| l. Again, it is difficult to reconcile this gap with the notion of a (half-) metallic ferromagnetic state below T c G3]. 
In half-metallic fcrromagnets, like CrC>2 or Fe3C>4, there is a band gap for electron states of only one spin direction. 
The opposite spin electrons have no gap at the Fermi level, similar to a standard metal situation. These states will 
contribute to tunnel current as in conventional metals, so that there would be no such temperature dependent gap 
feature in the tunnel spectroscopy H] like the one observed for the doped manganites [p3. 

We note that within the framework of our theory there should be a temperature dependent gap A related to the 
breakdown of a bipolaron into two polaronic carriers. The density of bipolarons peaks at T c , whereas the polaron 
density dips there (Figs. 1, 2) and, therefore, the gap feature in the tunneling I — V curves will be most pronounced 
in this region, as observed J24J ] . Spin-polarized polarons will provide a gapless background for tunneling current, 
which is least important in the vicinity of the transition temperature. We note that STM should also be sensitive 
to the presence of the one-particle charge-transfer gap between filled Mn d and empty O p states. In addition to the 
temperature dependence, we can predict how the gap feature will depend on the critical temperature of the transition 
(Fig. [}]). Namely, as already follows from our discussion, with the increase of A the critical temperature T c goes down 
[p8[ . Very similar behavior has indeed been observed experimentally on samples with different critical temperatures 
(Fig. §, inset) [Q. 

The giant isotope effect in Lao.sCao^MnOs, where a shift of -21K in T c was observed as a result of 16 to 18 
substitution pq ], is quantitatively explained within our approach. Namely, the gap is given by fl5| ] 

A = 2E P — Vc — ^W, (30) 

where E p is the polaron level shift, Vc is the Coulomb repulsion between bound polarons, and W — Wo exp(— g 2 ) is 
the polaron bandwidth rcnormalizcd from the bare value Wo with the electron-phonon interaction constant g 2 fig] . 
The only quantity in ( |30| ) that depends on ionic mass is the polaronic exponent g 2 — 7 E v / '(fiui) oc M 1 / 2 Ea], where 
7 < 1 is a numerical coefficient depending on the radius of the electron-phonon interaction |35| . As immediately 
follows from this relation, isotope substitution will change the gap A in the following way 

A 18 = A 16 + W^ 6 (Vl8/16-l), (31) 
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FIG. 7. Isotope effect on magnetization (a) and inverse carrier density (b) of Lao.8Cao.2Mn03+y calculated in the present 
theory. Inset: experimental results [16]. Substitution 16 0-^ 18 leads to increased resistivity (b). 



where indices mark the quantities for the corresponding isotopes of oxygen. According to ( |3l| ) Aig is always larger 
than A 16 . This automatically leads to a lowering of T c as a result of the isotope substitution, as observed |16|, in 
Fig. [r| The resistivity on the other hand, is larger in 18 substituted samples, and this correlation seems to be 
supported by recent experiments [ [f6| |. Note that the single parameter defining the isotope effect on the magnetic 
transition and the resistivity jump is Wg\§, since neither E p nor Vc depend on the ion mass. 



VIII. LOCALIZATION OF POLARONS BY DISORDER 

We have also studied the localization of p- holes due to a random field with a gap A/2 between localized impurity 
levels and the conduction band. The energy of polarons on impurity centers is given by 



E 



*TU) 



A 
~2~ 



1 



{+)-V ip m-(+)fx B H, 



(32) 



where Vi P is the exchange interaction between localized and delocalized polarons. The band diagram for this case is 
the same as in Fig. [ljwith the replacement of the bipolarons by localized polarons. 

Assuming that the Hubbard repulsion prevents a double occupancy of the impurity centers, one can easily obtain 
the thermodynamic potential for impurities 



-fcsTln 



1 + 2uye A / 2kBT cosh 



W ip m + /j, b H\ 



The chemical potential is found to be 



2un cosh(C) 



exp(-S/t), 



(33) 



(34) 



where £ = (^Vi p m + HbH) /k B T , if we assume that the total number of impurity states is x{= doping). 

We have found similar features of the phase transition in zero field in the impurity case, as compared with the 
previous case with bipolarons. Thus, we obtain, by linearizing the system of equations of state (|l6|)-(|l§|) with (|34"|), 
the following equation for the polaron density at the transition in zero field: 



^ln: 



= 2 1 ' 2 5. 



(35) 
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Apparently, it has solutions only for S < S c (x). Therefore, the transition is first order for 5 = A/J p dS > 5 c (x) and 
second order for <5 < S c (x), with S c (x) slightly larger than in the case of the bipolaron localization. This follows from 
the same consideration as in our previous discussion of Eq. ( |2l| ) . 

The field sensitivity in the case of disorder localized polarons is much lower than for the bipolarons. This stems 
from the different functional dependence of the chemical potential. The present approximations are valid in the limit 
y -C 1, meaning that the polaron carriers are non-degenerate. In contrast to the case of the bound polaron pair 
formation, in the impurity case the expression ([54j) for y is singular, y cx 1/n, in the limit of small polaron density. 
This means that in the vicinity of the current carrier density collapse the value of y sharply increases in the case of 
polarons localized on impurities. As a result, the collapse becomes less pronounced, and transport becomes far less 
sensitive to an external field. We note also that Eq. ( p4|) contains a factor depending on the external magnetic field 
in the denominator. This is in contrast with the case of bipolarons (^8|), where the field dependence is suppressed 
by a small factor exp(— J st //csT c ). This field dependence, however, is small since always ixbH I^bTc "C 1, and it 
quickly vanishes in the low-temperature phase when the exchange interaction sets in (m ^ 0), as one can see from 
the expression for the parameter £ above. The singular behavior of y as a function of the density for n — ► 0, and the 
Zeeman splitting of the impurity states makes the transition far less sensitive to the magnetic field. As a result, no 
quantitative description of the experimental CMR data has been found with the localization of polarons due to disorder. 



IX. CONCLUSION 



In conclusion, we have developed a theory of the ferromagnetic-paramagnetic phase transition in doped magnetic 
charge-transfer insulators with a strong electron-phonon coupling. We have found that a few non-degenerate polarons 
in the p band polarize localized d electrons because of the huge density of states in the narrow polaronic band. For 
a sufficiently large p — d exchange J pe iS > A, we have obtained a current carrier density collapse at the transition 
owing to the formation of immobile local pairs in the paramagnetic phase with the binding energy A about twice that 
of the polaron level shift (IB). Depending on the ratio A/(J p dS), the transition is first or second order fl38f . 

We have explained the resistivity peak and the colossal magnetoresistance of doped perovskite manganites, Fig. 0, 
as the result of the current carrier density collapse due to the binding of polarons into local pairs (bipolarons). The 
density of these immobile pairs has a sharp peak at the ferromagnetic transition when the system is cooled down 
through the critical temperature T c . Below T c the binding of polarons into pairs competes with the ferromagnetic 
exchange of p-holes with the Mn d A local moments, which tends to align the polaron moments and, therefore, breaks 
those pairs apart. The spin-polarized polaron band falls below the bipolaron band upon decrease in temperature, so 
that all carriers are unpaired at T = if J p dS > A. Above T c , the bipolaron density decreases because of thermal 
activation across the polaron binding energy. These competing interactions lead to the unusual behavior of CMR 
materials, the huge sensitivity of their transport to external field, and the very large negative magnetoresistance. 

There is a crossover around the transition temperature from polaron tunneling at low temperatures to polaron 
hopping, where the latter dominates at high temperatures. This explains the temperature behavior of the resistivity 
in a wide temperature range around the transition. The ferromagnetic to paramagnetic transition is also accompanied 
by a sharp anomaly in the specific heat. 

The present theory provides a natural explanation for the temperature dependent gap feature in tunneling spectra 
[ p4| and the giant isotope effect on the temperature of the ferromagnetic transition . One of our main conclusions 
is that the highly polarized ferromagnetic phase of manganites is a polaronic doped semiconductor rather than a 
metal. 

We expect that the present theory is general enough to also account for the giant magnetoresistance observed 
in pyrochlore manganites (4^] and other systems ||. It is worth mentioning in this regard that the present theory 
requires the presence of strong electron-phonon coupling of any origin, but it does not require the presence of Jahn- 
Teller distortions and/or the double exchange mechanism. Note that the Jahn- Teller distortions and the double 
exchange mechanism are certainly absent in, for instance, pyrochlore manganites, chromium spinels and other 
CMR systems, so that the ideas based on the double exchange cannot be applied there at all. It is believed that at 
least in perovskite manganites the local Jahn- Teller distortion may be involved in defining the crystal structure of the 
parent insulating phases [Q, although tilting distortions of MnOg octahedra are just a result of steric conditions [ff9|||. 
Apparently, the ratio of the sum of Mn and O ionic radii, rMn+J'o, and (rL a + ?'o)/v / 2 (misfit parameter) substantially 
differs from unity to make a cubic structure unstable and favor a rotation of MnC-6 octahedra pp9|] . The tetragonal 
distortion of MnOg is large, its symmetry corresponds to a notion of the Jahn- Teller local distortion. However, since 
the steric interaction is strong, it necessarily deforms the lattice, thus rendering the Jahn- Teller derivation, strictly 
speaking, inapplicable. 

It is also believed that doping by divalent metals introduces holes into the Mn 3+ d-shell, since the doped systems 
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are less distorted |50). This argument, which may have supported the relevance of the double exchange mechanism for 
at least perovskite manganites, contradicts the site-sensitive spectroscopic probes |p2|p6| , which show unambiguously 
that holes reside on O sites. It also neglects two important facts, that (i) the doping is heavy (> 10 21 e/cm 3 ) and 
there is a substantial size difference between the impurity and host atoms and (ii) the O p-holes are hybridized 
with the d states on Mn 3+ , depending on the value of the charge-transfer gap. Both effects, together with screened 
Coulomb holc-holc repulsion, can apparently explain the observed changes in the lattice distortion upon doping 
without invoking the Jahn- Teller mechanism. These short-range interactions may well be responsible for the charge- 
ordered phases observed at some doping levels in manganites |9(]. It would be interesting, in this regard, to perform 
quantum-chemical calculations of MnC>6 clusters with holes doped onto O site(s). 

Changes and the amount of disorder in the bond lengths are very important for characterizing the properties 
of polaronic systems. The reduction in bond length distribution width as a result of cooling through T c in doped 
manganites has been attributed to (at least partial) delocalization of doped carriers in low-temperature 'metallic' 
phase. Since the data shows that the carriers retain their polaronic character below T Cl and the residual width of 
the Mn-0 bond length distribution remains larger than that of CaMn03 [p0[ , where the Jahn- Teller Mn 3+ ions are 
absent, the reduction of the width should be mainly related to instability of bipolarons in this temperature region. 
Breaking of polaron bound pairs below T c may result in a reduction of bond length distribution width, and we shall 
address this question elsewhere. It is worth repeating that whether or not the Jahn- Teller distortions play any role 
in doped perovskite manganites and the exact location of the carriers is of no importance for the present scenario of 
the CMR. 

We acknowledge useful discussions with A.R. Bishop, D.M. Edwards, J. P. Franck, K.M. Krishnan, P.B. Littlcwood, 
S. von Molnar, V.G. Orlov, W.E. Pickett, D.J. Singh, S.A. Trugman, and R.S. Williams. We especially grateful 
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